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Results are presented from numerical simulations of the flat-space nonlinear Klein-Gordon equa- 
tion with an asymmetric double-well potential in spherical symmetry. Exit criteria are denned for 
the simulations that are used to help understand the boundaries of the basins of attraction for 
Gaussian "bubble" initial data. The first exit criteria, based on the immediate collapse or expan- 
sion of bubble radius, is used to observe the departure of the scalar field from a static intermediate 
attractor solution. The boundary separating these two behaviors in parameter space is smooth and 
demonstrates a time-scaling law with an exponent that depends on the asymmetry of the potential. 
The second exit criteria differentiates between the creation of an expanding true- vacuum bubble and 
dispersion of the field leaving the false vacuum; the boundary separating these basins of attraction 
is shown to demonstrate fractal behavior. The basins are defined by the number of bounces that the 
field undergoes before inducing a phase transition. A third, hybrid exit criteria is used to determine 
the location of the boundary to arbitrary precision and to characterize the threshold behavior. The 
possible effects this behavior might have on cosmological phase transitions are briefly discussed. 
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I. INTRODUCTION 

While expanding and collapsing scalar field bub- 
bles have each been given a great deal of attention 
in the physics community (largely in the contexts of 
phase transitions and oscillons, respectively), there 
is surprisingly little focus on the solutions that live 
on the threshold of those two end states. This pa- 
per focuses on the critical phenomena that exist near 
this threshold and the fractal boundaries that define 
it. 

The model used in this paper is the nonlinear 
Klein-Gordon (nlKG) equation with an asymmet- 
ric double- well potential (ADWP). The ADWP used 
here is the typical quartic potential with two stable 
minima: one global minimum at the lowest possible 
potential energy (the "true" vacuum) and one local 
minimum at a higher energy (the "false" vacuum). 
In the study of phase transitions, one typically con- 
siders initial data that consist of perturbations on 
top of the false vacuum. If the perturbations are 
large enough, the tendency for the field to transition 
to the energetically more favorable true vacuum can 
overcome the effective "surface tension" that would 
cause the bubble to collapse. When that is the case, 
the bubble wall will expand outward causing a tran- 
sition from the false to the true vacuum; the result- 
ing solution is referred to as a supercritical bubble. If 
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the perturbations are small, however, a phase transi- 
tion will not occur; the scalar field will collapse and 
eventually disperse, and the solution is referred to 
as a subcritical bubble. 

When the nlKG equation was first studied in the 
context of phase transitions, research efforts focused 
largely on supercritical bubbles, and in particular, 
on domain walls. Subcritical bubbles tended to be 
neglected because, by definition, they do not induce 
a phase transition; it was not understood until much 
later that subcritical bubbles could still have cosmo- 
logical implications. For a wide range of initial data 
leading to subcritical bubbles, there exist localized 
oscillating solutions, called oscillons, that have life- 
times that are large compared to naive expectations. 
Oscillons were originally discovered by Bogolubsky 
and Makhankov [1] and have been thoroughly in- 
vestigated by numerous authors using many types 
of nonlinear models [2rTT2] ; the effect oscillons may 
have on phase transitions has been of great interest 
[13HT5] . The work by Gleiser and Sicilia [16] is an 
excellent summary and analytic characterization of 
the salient features of oscillons. 

While most of the simulations presented in this 
paper result in the formation of either an oscillon or 
a domain wall, this paper actually focuses on nei- 
ther; this paper focuses on the nonlinear dynamics 
that exist on the threshold of creating one of these 
two end states. This paper begins with a presenta- 
tion of general formalism and definitions in Section 
|H| In Section |III| methods are presented that can 
be used to solve for the intermediate attractor solu- 
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tions to the nlKG-ADWP system. It is shown that 
in the neighborhood of the attractor, Type I crit- 
ical phenomena similar to those observed in grav- 
itational collapse are present. Section [TV] demon- 
strates that the basin boundaries have fractal struc- 
ture and can be parameterized by the number of 
times they "bounce" before inducing a phase tran- 
sition, not unlike the kink/antikink soliton (domain 
wall) collisions of ]T7]. Furthermore, methods are 
presented that allow for the arbitrarily precise de- 
termination of the location of the fractal edges. Be- 
ing able to explore these edges in parameter space 
reveals that both a static intermediate attractor so- 
lution and Type I critical phenomena appear to ex- 
ist on every basin boundary. Finally, Appendix [A] 
presents a rough approximation for the Minkowski- 
Bouligand fractal dimension of a data set with the 
basic properties observed in the nlKG-ADWP sys- 
tem. 



II. GENERAL FORMALISM AND 
DEFINITIONS 



conserved energy 



The nlKG action is 



S[<t>] = / dWM --sTd^W-vty) , (i) 



where <j) = (j)(r, t) and g^ v is the flat spacetime met- 
ric on a four-dimensional manifold. The potential V 
is given by V{4>) = \<t> 2 ((j) 2 - 4(1 + 8)<j) + 4), where 
6 is a measure of the asymmetry of the potential. 
When (5 = 0, V((f>) is a symmetric double-well po- 
tential with vacuum states at = and <f> = 2. With 
5 7^ 0, the false and true vacuums are fa = and 
fa = |(1 + 8) + |Vl + 18(5 + S 2 , respectively. 

Variation of the action with respect to cj> gives the 
covariant equation of motion 

-^-d, U\g\9^d v A = <t> (0 2 - 3(1 + (5)0 + 2) , 

(2) 

which after imposing spherical symmetry and using 
spherical coordinates (i,r) gives 

n = ^d r {r 2 t>) -0(0 2 -3(l + (5)0 + 2) (3) 



4> = n, 



(4) 
(5) 



where II = dt4> and $ r = d r (j>. The spacetime ad- 
mits a timelike Killing vector that gives rise to the 



E = Air 



dr 'i 



(n 2 + $ 2 ) + v(<f>) - v(fa) 



where is the radius of the outer boundary, and 
the constant V(fa) is subtracted to give E = for 
(f>(r,t) — fa (i-e., everywhere at the true vacuum). 

The initial data used to evolve ([3|, (4J), and Q 
are taken to be a simple Gaussian bubble: 



4>{r, 0) = fa + {fa - <M exp (-r 2 /cr 2 ) 



(7) 



This bubble type of initial data interpolates between 
the true and false vacuums and represents a pertur- 
bation away from the false vacuum that may (or may 
not) lead to a phase transition. 

To explore the dynamics of the scalar field bub- 
bles, it is useful to define an ancillary function 



m) = 



max(r, (fa F ,tj)) when 0tf € <H r i> tj) 
otherwise, 

where fa F = (fa + fa)/2, and fap satisfies fa F £ 
(f){n,tj) at tj if &{n,tj) < fa F < <j>{r i+ i,tj) for 
some i. More simply put, £(t) is the maximum ra- 
dius for which the field is halfway between the true 
and false vacuums, and zero if at time tj the field 



us- 



does not anywhere equal fa F . For example, 
ing the initial data profile ([7| and the definition (|8| , 
one can define £o = £{to = 0), which is simply the 
radius (at half-maximum) of the Gaussian profile, 
£o = croVln 2. For supercritical bubbles, £(i) is a rea- 
sonable approximation to the location of the bubble 
wall. 

For collapsing bubbles, the number of oscillations, 
or bounces, of £(t), denoted N (£ (t j), is defined to 
be the number of times 



> 



and 
0, 



(9) 



for any j. To retain common terminology with the 
works on kink/antikink soliton collisions, the term 
"bounce" will be used to describe the oscillations 
the field undergoes. 

Finally, all solutions to ([3]), Q, and ^ are 
obtained using a second-order iterative Crank- 
Nicholson finite difference scheme with higher-order 
Kreiss-Oliger dissipation [2J. A simple outgoing 
boundary condition is used at the exterior of the 
grid, but the grid is chosen to be large enough such 
that any reflected radiation could not travel back to 
origin during the evolution. Since this work is fo- 
cused on the threshold of bubble expansion and not 
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FIG. 1: Bitmap of two-dimensional parameter space 
survey showing the collapse (black region) or expansion 
(white region) of initial data with Gaussian width oo 
and asymmetry parameter S. £ i is used to determine 
expansion or collapse for 62 500 (250x250) pairs (ao,8), 
for 1.0 < (To < 7.5 and 0.05 < 6 < 0.25. The boundary 
between collapse and expansion is smooth. 



the long-term behavior of oscillons, this approach is 
simple and effective. 
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FIG. 2: Plot of solution lifetime as a function of 
— In |oo — 0ol for 8 — 0.141 using exit criteria £%. Time- 
scaling exponent 7 = 1.62 is observed for both expanding 
(top line) and collapsing (bottom line) evolutions. Sim- 
ilar scaling exponents 7 = {3.94,2.13,1.54,1.23,1.04} 
were measured for 5 = {0.05,0.10,0.15,0.20,0.25}, re- 
spectively. 



A. Time-scaling law 



III. THE INTERMEDIATE ATTRACTOR 

It has been shown in [IB] that solutions to ([3| , Q , 
and ([5| obey a time-scaling law near the threshold of 
expansion or collapse. The exit criteria for the sim- 
ulations in [18| is based on the bubble radius moving 
sufficiently far away from the initial bubble radius, 
£0: 



£1 



expansion when £(t) > £ + 
collapse when £(t) < £~ , 



(10) 



where £ + and £ - are constants chosen to be above 
and below £o> respectively. Figure [T] displays the end 
state (expansion or collapse) as a function of the 
width of the initial bubble and the asymmetry 
parameter 5. The boundary between the two end 
states is observed to be smooth, and the <5-dependent 
threshold value of Co is defined to be <Tq(S), or just 



To explore the properties of the threshold, one 
must find <7q precisely. This is done by fixing the 
asymmetry parameter S and varying <7o- The search 
begins with a Og~ that is known to expand and a <7q 
that is known to collapse. Evolutions are run with 
test values that successively bisect <7q~ and <Tq until 
machine precision is reached (64-bit double-precision 
variables used). For a fixed asymmetry parameter 
5, one can observe the lifetime of the solution as a 
function of a®. The solution lifetimes are observed 
to obey a time-scaling law 



T = -7ln|cr - <Tq 



(11) 



with an exponent 7 that depends on <5. The time- 
scaling law can be observed on both the expansion 
and collapse sides of o-q in Fig. ]2| 

The time scaling laws like (11) are not uncom- 



mon in the study of nonlinear field theories. Similar 
Type I critical phenomena are observed in the grav- 
itational collapse of a wide range of matter fields 
[T§H25]. In these models, in the limit that ao ap- 
proaches <7q, the solutions each approach an unsta- 
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FIG. 3: (Color figure and movie online) Time evolu- 
tion of 4>{r,t) for static (solid black), collapsing (dot- 
ted green), and expanding (dashed red) solutions with 
S = 0.141, using E\. Initial bubble radius oq is fine-tuned 
to o"q to within machine precision (64-bit floating point). 
The collapsing and expanding solutions oscillate about 
the static solution and disperse, leaving the spacetime in 
different vacuum states. 



ble intermediate attractor solution with Lyapunov 
exponent I/7. While the exact nature of the inter- 
mediate attractor solutions varies depending on the 
physical model, the most common attractor type is 
static [5TH24] . 



B. The nlKG static intermediate attractor 

Not surprisingly, observation of the threshold dy- 
namic solutions shows that the field appears to con- 
sist of small oscillations on top of an underlying 
static solution. Fig. [3] shows snapshots of a near- 
critical evolution for values of <7o slightly above and 
below (Tq . Guided by this observation and other ex- 
amples of Type I critical phenomena, a static solu- 
tion to ([3]), Q, and ([5]) is obtained by invoking a 
static ansatz (setting all time derivatives equal to 
zero). The static field equations are 

\d r (r 2 $) = 0(0 2 -3(l + £)0 + 2) (12) 



FIG. 4: (Color online) Solution <f>*(r) obtained from the 
solution of Q and Q with S = 0.141. This solution is 
the ^-dependent static intermediate attractor for (4b, 
and |5| at the collapse/expansion threshold fine-tuned to 
machine precision. The "x" denotes the value of £ w 2.3. 
The inset shows —V((f>). 



d r (f> = $, 



(13) 



and are solved by "shooting." 

The solution of (fl2l and ( 13 ) is very similar to 



the one-dimensional particle subject to a potential 
of — V{4>), differing only by the geometric factor in 
the spatial Laplacian and the exchange of r and t. 



The critical solution is obtained by solving ( 12 ) and 
(13) with 3>(r = 0) = and different initial values 
of (j>(r — 0). The solutions have three types of be- 
havior as r — > 00: 4> —¥ ±00 or oscillation in the 
local minimum of — V((f>). Bisecting on the thresh- 
old of the <f> — > —00 end state and the oscillation end 
state yields the desired critical solution that asymp- 
totically approaches the false vacuum. In the one- 
dimensional particle analogy, this is like finding the 
perfect initial position of the particle so it can "roll" 
through the trough of the potential and come to a 
perfect stop at the unstable minimum at (f> = (Fig. 
[4| inset). Fig. BJshows the fine-tuned solution to ( 12 ) 
and (pi for 5 = 0.141. 



In addition to showing the evolution of the dy- 
namic solutions, Fig. [3] also shows the static in- 
termediate attractor solution. It appears that the 
threshold solution is a superposition of the static so- 
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FIG. 5: (Color online) Plot of obtained from evolv- 
ing Q, and |5| for two values of <to, one slightly 
above oq and one slightly below for S = 0.141. The so- 
lution with (To > gq leads to an expanding true-vacuum 
bubble. The solution with ao < <tq leads to a bubble 
that collapses. The horizontal line represents the value 
of £ for the static intermediate attractor. 



FIG. 6: (Color online) Plot of £(t) for three solutions 
subject to criteria £2- Solution A is supercritical. Solu- 
tion B is subcritical and exits due to the N(£(t)) > iVosc 
condition (only three oscillations are shown). Solution 
C is subcritical and exits due to the t > t m ax condition. 
Since solution C has so little energy, it does not oscillate 
even once, i.e., N(£(t)) = for all t. 



lution and a small amplitude time-dependent local- 
ized "shape mode," reminiscent of one-dimensional 
kink/antikink soliton collisions [25]. Fig. [5] shows 
the approximate bubble radius, £(t), for the same 
two near-critical solutions shown in Fig. |3j it can 
be clearly seen that the bubble radius for the dy- 
namic solutions oscillates about the bubble radius 
of the static intermediate attractor before collapsing 
or expanding. 



IV. TRUE/FALSE VACUUM BOUNDARIES 

In the previous section it was shown that E\ di- 
vides the ag-S parameter space into regions of imme- 
diate expansion or collapse. This section focuses on 
the more physical concern of what vacuum state the 
spacetime approaches in the t — > 00 limit. In regions 
where the bubbles expand, £ > £ + , the bubble wall 
expands outward indefinitely and induces a phase 
transition to the true vacuum. In regions where the 
bubbles collapse, however, E\ terminates the evolu- 
tion immediately after drops below £ - , and as 
such, no conclusion can be drawn about which vac- 



uum state the field will occupy as t — > 00. There are 
three long-term behaviors observed in the collapsing 
field configurations: solutions that immediately dis- 
perse to leave the false vacuum, solutions that form 
oscillons that eventualy disperse to leave the false 
vacuum, and solutions that collapse but bounce back 
and lead to an expanding true-vacuum bubble. 

An exit criteria that effectively determines which 
vacuum state the solution will approach as t — > 00 
is 

!PT when £(t) > 
NPT when JV(£(t)) > N osc (14) 
or t > i max , 

where PT and NPT refer to phase transition and 
no phase transition, respectively; £ + is again the ra- 
dius beyond which the bubble is assumed to induce a 
phase transition; A^ osc is the number of bounces be- 
fore assuming the solution enters the oscillon state; 
and i max is a time at which the bubble is assumed 
not to induce a phase transition (if £(t) has remained 
less than 

Figure[6]displays three solutions subject to criteria 
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FIG. 7: (Color online) Bitmaps of two-dimensional pa- 
rameter space survey showing the collapse (black region) 
or expansion (white or colored regions) of initial data 
with Gaussian width ao (horizontal axis) and asymme- 
try parameter 5 (vertical axis). £2 is used to determine 
expansion/collapse, and different shades of grey (colors 
online) denote PT solutions with different numbers of 
bounces. The top bitmap displays 1.0 < 00 < 7.5 and 
0.05 < S < 0.25, and the boxes in each bitmap denote 
the field of view of the next bitmap. 
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FIG. 8: (Color online) Plot of {(t) for two solutions 
subject to criteria £2- Solution A is PT but exits with 
N (£(*)) = 3 (i- e -> three bounces). Solution B is NPT 
and exits due to the iV(£(t)) > N OBC condition. 



£2, each exiting due to a different condition (denoted 
A, B, and C in the figure). Solution A is an exam- 
ple of a typical runaway PT solution; these solutions 
are easily detected with the £(t) > £ + condition and 
a reasonably chosen ^-dependent £ + . Solution B is 
an example of a solution that settles into the oscil- 
lon regime, staying localized for many periods before 
eventually dispersing (NPT); these solutions are de- 
tected with the N(£(t)) > N osc condition. The vast 
majority of the solutions discussed in this paper are 
of these two types (A or B). The t > i max condition 
is added to the NPT criteria for the cases like solu- 
tion C, where the initial data perturbations are so 
small that the field immediately disperses without 
enough energy to even form an oscillon. 

Figure [7| shows a bitmap displaying the PT/NPT 
end state as a function of a and S for solutions 
subject to £2- A remarkable difference can be seen 
when compared to Fig. [I] The simple boundary be- 
tween collapsing and expanding solutions is replaced 
by a much more complicated fractal boundary. The 
extra structure is a result of solutions that bounce 
back after collapse to form expanding true-vacuum 
bubbles (see Fig. [8]). These solutions collapse ini- 
tially yet still lead to phase transitions. The solu- 
tions shown in Fig. 7] are shaded (colored online) dif- 
ferently based on the number of times they bounce 
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before expanding to induce a phase transition. Each 
band of n bounces is immediately surrounded by a 
set of bands of (n+1) bounces that approach it from 
both sides in (j -space. The structure appears to be 
fractal given it continues to repeat itself in a self- 
similar fashion upon successive magnification. This 
structure closely resembles the structure observed 
in [T7] , except that instead of colliding domain walls 
(kink/antikink soliton collisions) parameterized by 
boost velocity, it arises from collapsing bubbles pa- 
rameterized by initial bubble radius. 

The width of er -space surrounding Cg that 
demonstrates this fractal bounce behavior is defined 
to be Acrg; given this, the region of <7o-space demon- 
strating this behavior is given by 

S* = {cr : (oq - A<7o) < (T < Oq}. (15) 

In Fig. [7J it is observed that for S « 0.141, CTq ps 2.52 
and Ao-q « 0.07. 



A. Threshold of the re-bounce band 




FIG. 9: Plot of solution lifetime as a function of 
— ln|o"o — <5"o| for ao on the PT side of <tq subject to 
exit criteria £2, each with different <5o values obtained 
by using £3 with different iV m i n , for asymmetry param- 
eter 8 — 0.141. The time-scaling law is observed until 
the ao values become large enough to lie outside the PT 
band being observed (structure on left side of plots). 



After observing the fractal structure of the 
PT/NPT boundaries in Fig. [7j one might wonder 
whether there exist scaling laws at the threshold of 
each boundary, and if so, whether the scaling law 
exponents and intermediate attractor solutions are 
universal. To answer this question, one must be able 
to determine the boundary of a given n-bounce band 
to arbitrary precision. One might attempt to use 



the same bisecting technique used in Section III A 
but with the PT/NPT condition £ 2 . Unfortunately, 
this technique fails due to the fractal nature of the 
boundary. For a fixed 5, one cannot use £2 to find 
a a PT (PT end state) and a a^ PT (NPT end state) 
that bound a single PT/NPT boundary. In fact, it 
appears that there would be an infinite number of 
such boundaries. Any bisecting approach using £2 
that began with a a PT in the desired PT band and a 
(Jq PT outside the band would almost definitely lead 
to a PT hopping to a different band. 

However, one can find the boundary of an n- 
bounce band if one exploits the fact that an n- 
bounce band is approached by either NPT solutions 
that form oscillons or PT solutions with more than n 
bounces. If one starts with a a PT that is in the A^ m i n - 
bounce band and a a® that is sufficiently close but 



not in the band, the hybrid condition 



PT 



when £(t) > 



^ delayed-collapse when N(£(t)) > A min 
and £(*) < £" 

(16) 

will successfully find the desired edge, denoted Oq, 
where ctq 6 S* Solutions inside the 7V m i n -bounce 
band will induce a phase transition and exit due to 
the £(f) > £ + condition. Solutions outside the -/V min - 
bounce band will bounce more than N m - m times since 
they either form an oscillon or are contained in a PT 
band with more than iV m i n bounces; when these so- 
lutions bounce more than N m - U1 times, they exit as 
what is here referred to as a delayed-collapse solu- 
tion. 

To understand the nature of solutions near the 
fractal boundary, solutions are first evolved with exit 
criteria £ 3 in order to find cJq . Then, solutions are 
obtained with exit criteria £ 2 using ao that approach 
Ctq from either side, so that one can observe the PT 
or NPT end state of the solutions near an 
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FIG. 10: Plot of solution lifetime as a function of 
— ln|ero — oq| for <ro on the delayed-collapse side of 
subject to exit criteria £2, each with different erg values 
obtained by using £3 with different AT m i n , for asymme- 
try parameter S = 0.141. Each iV m i n threshold solution 
shows roughly equal-width bands of PT solutions ap- 
proaching ctq in roughly log-periodic windows of param- 
eter space. 



FIG. 11: Plot of £(t) for PT and delayed-collapse so- 
lutions for four Af m i n values, S = 0.141, and 00 values 
immediately above and below four distinct critical solu- 
tions (one for each 7V m i n ). The horizontal lines in each 
plot correspond to the £ value for the S = 0.141 static 
intermediate attractor. Animations of numerical evolu- 
tions are available online [? 1. 



demonstrates the time-scaling law 
T = -7ln|<7 - &o\ 



(17) 



for (To approaching (Tq on the PT side. While similar 
to the time-scaling seen in Section [ill A| there now 
appears to be a scaling law observed on each of the 
infinitely many PT/NPT boundaries. Furthermore, 
the exponent 7 appears to be universal for a given 
5. Fig. [T0| shows a mix of PT and NPT solutions for 
<7o approaching <7q from the delayed-collapse side. 
The figure seems to imply that each band of PT so- 
lutions with n bounces is approached by an infinite 
number of PT bands of (n + 1) bounces. For ex- 
ample, in the N m [ n = 1 plot in Fig. [Toj each of 
the curves represent PT solutions with N(£(t)) > 2 
bounces. The main dips that can be seen are for PT 
solutions with N(£(t)) = 2 bounces; the regions of 
fine structure seen around the N(£(t)) = 2 bands 
each have iV (£(£)) > 2 bounces, showing that the 
structure continues at finer and finer scales. The 
gaps between the bands are for 00 that form oscil- 
lons and eventually disperse (with no phase tran- 
sition). While the appearance of the bands might 



only be a truncated fractal structure, the fact that 
the bands appear in roughly log-periodic windows 
that approach <7q to machine precision highly sug- 
gests that the bands would continue ad infinitum 
(Appendix [A] discusses approximations to the frac- 
tal dimension). The time-scaling law is present yet 
again (with the same (5-dependent 7) but here serves 
as a lower-limit envelope on the lifetime. 

(animations online, 



Finally, looking at Fig. 
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]) one can see that solutions on this fractal boundary 
indeed approach the (S-dependent static intermediate 
attractor. This strongly suggests that the solution 
obtained from the static ansatz (for a given 5) serves 
as a universal intermediate attractor on the thresh- 
old of inducing a phase transition. While the frac- 
tal structure itself is not surprising given the previ- 
ous understanding of kink/antikink soliton collisions 
[TT] . observing this structure in bubble collapse, the 
direct measuring of the time scaling law, the ability 
to determine the fractal boundary <7q to arbitrary 
precision, and the realization that the static solution 



to (12 1 and (131 serves as an intermediate attractor 



on the fractal boundary all offer new insight into a 
classic nonlinear system. 
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V. DISCUSSION AND CONCLUSIONS 

It has been shown that the boundary between 
the basins of attraction for the nlKG-ADWP sys- 
tem (the true- and false-vacuum states) is fractal in 
nature. The fractal structure arises from solutions 
that collapse, bounce a number of times, and then 
expand to induce a phase transition. The solutions 
that induce a phase transition are organized into 
bands that are grouped by the number of bounces 
they undergo before becoming expanding bubbles. 
These bands form a regular and predictable (as op- 
posed to chaotic) fractal structure in the a — 6 pa- 
rameter space. Bands of solutions with n bounces 
are surrounded by bands of solutions with (n + 1) 
bounces. When viewed in ln|<7o — <Jo|-space, the 
roughly equal-width (n + l)-bounce bands appear 
to approach the edges of the n-bounce bands in a 
log-periodic fashion. A time-scaling law is observed 
at each fractal boundary with an exponent that de- 
pends on the asymmetry parameter <5, and the so- 
lutions that exist on these boundaries are shown to 
approach the same static intermediate attractor so- 
lution. 

While the phenomena discussed here have been 
studied in the context of isolated scalar field col- 
lapse, one may wonder what possible impact these 
phenomena would have on global phase transitions 
induced by bubble nucleation. The Gaussian bub- 
bles Q are assumed to represent perturbations away 
from a false vacuum. For a global phase transition 
to occur, supercritical perturbations need to form 
expanding bubbles that combine, thereby inducing 
a global transition to the true vacuum. The impact 
the fractal behavior discussed in this paper might ac- 
tually have on bubble nucleation would depend not 
only on how close to threshold these perturbations 
are for a given system, but also on how those per- 
turbations vary throughout space. For example, if 
typical perturbations are near the threshold of ex- 
pansion/collapse, (<7o) ~ (Jq, but the variation in ao 
in different perturbations throughout space is large, 
Ado > AiTq, the probability that these effects will 
have an impact on previously explored bubble nu- 
cleation rates will be low since most ao are not con- 
tained in E* and therefore do not demonstrate the 
fractal bounce behavior. However, if Aao < Actq, 
then most ao are in £*, and the phenomena ob- 
served in this paper would likely have a significant 
impact on bubble nucleation rates. In this condition, 
one might want a way to approximate the probabil- 
ity that a given perturbation (given by ao) would 
end up in either the true- or false- vacuum state; one 
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FIG. 12: Plots demonstrating two observed fractal prop- 
erties of the PT solutions. When viewed in In \ ao — <Tq|- 
space, the PT solutions appear in roughly equal-width 
(n + 1) -bounce bands that approach the edges of the 
n-bounce bands in a log-periodic fashion. The top and 
bottom plots show this behavior as a function of ao and 
— hi | ctq — o"ol) respectively, for oq — 2.5. 



very natural measure of this probability is related 
to the fractal dimension itself, which is discussed 
in Appendix A. Furthermore, one would not need 
a stochastic description of perturbations to create 
a phase transition that appears chaotic or random; 
a smoothly varying distribution of perturbations, 
where ao £ £*, would still likely create a system 
with a rich and complex structure. 



Appendix A: Minkowski-Bouligand dimension of 
log-periodic line segments 

The structure observed in Fig. [7] seems to be 
fractal since a self-similar banded structure appears 
at successive scales of magnification. For each n- 
bounce band, there appear to be an infinite number 
of (n + l)-bounce bands of roughly equal width in 
In | (To ~ <To I -space that approach the n-bounce band 
on either side of cro-space in a log-periodic fashion. 
By defining a set on R 1 that demonstrates similar be- 
haviors, this appendix calculates an approximation 
to the fractal dimension of the PT/NPT boundary 
discussed in this paper. 
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FIG. 13: Plots of the Minkowski-Bouligand dimension 
of log-periodic line segments for different values of A2 as 
a function of A1/A2. 



The notion of dimensionality used here is the 
Minkowski-Bouligand, or box-counting, dimension, 



D 



M 



lim 

e->0 



logJV(e) 

-log 6 



(Al) 



where N(e) is the number of boxes of side length e 
required to cover the set of PT solutions. In calcu- 
lating Dm, one can start with a single box covering 
the region in question. One then successively breaks 
the box into smaller boxes. The boxes that cover 
the set contribute to N(e), and the dimension Dm 
is obtained in the limit e — > 0. After fixing 6, the 
ao-5 boxes reduce to line segments parameterized by 
<7q. Choosing a box that is in the bulk of a PT band 



yields Dm — 1, since every box contains PT solu- 
tions and is in the set, no matter how many times it 
is divided into smaller boxes. Choosing a box that 
is far away from a PT band and in a clearly NPT 
region has Dm = 0, since no smaller box ever con- 
tains a PT solution allowing it to belong to the set. 
The intermediate region is a different story entirely. 



Figure 12 demonstrates the two observed proper- 
ties of the PT solutions. The first property is that 
the bands of PT solutions appear to be roughly equal 
width in ln|<7o — CTpl-space, here denoted Ai. The 
second property is that the bands of PT solutions ap- 
pear to be an equal distance, A2, apart in In |cr — 6-q |- 
space. The ratio Ai/A 2 is essentially a duty cycle 
quantifying what proportion of the <7rj space is in the 
set of PT solutions. It is this observation that would 
allow one with limited measurements of Ai and A2 
to approximate the fractal dimension of the system; 
these numbers could be used to determine the prob- 
ability that any <ro G X* will induce a phase tran- 
sition. Figure [13] shows the Minkowski-Bouligand 
dimension for different values of A 2 as a function of 
Ai/A 2 . As expected, all solutions have Dm = for 
Ai/A 2 = since those sets are empty. Likewise, all 
solutions have Dm = 1 for A1/A2 = 1 since those 
sets are completely solid lines. 

The approximation discussed here is only a very 
rough approximation to the fractal dimension ob- 
served in the nlKG-ADWP dynamical system since 
it only considers a single level of recursion of the 
self-similar pattern. In the actual nlKG-ADWP sys- 
tem, in addition to each n-bounce window being ap- 
proached by an infinite number of (n + l)-bounce 
windows, each (n + l)~bounce window also is ap- 
proached by an infinite number of (n + 2)-bounce 
windows, and so on. Therefore, for an observed Ai 
and A 2 , the values reported here should be consid- 
ered a lower bound. 
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